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Abstract 

A system of 'telegrapher's' equations for a number of long lines joined into a 
network is studied. Explicit formulas for Fourier transforms of current and voltage 
are derived. These formulas are very suitable for computer application as well as 
for the analytical study of processes o networks. As an example, the availability of 
formulas aids the derivation of explicit formulas for maxima of current amplitude 
over the given class of admissible external influences. These values may be used to 
indicate the characteristic of network robustness to excess voltage or electromag- 
netic impulse. The approach is based on the operational solution already proposed 
by the author for more general partial differential equations on graphs. 
Indexing terms: Telegrapher's equations, Branching long lines, Waves on graphs 

1 Introduction 

The classical 'telegrapher's' equation describes the evolution of current in a long line 
[1]. There are a lot of technical and physical objects that are modeled by a system (net- 
work) of jointed long lines( a network of underground pipes, a system of grounding, an 
ordinary electrical network of two-wire transmission lines with consumers and sources 
etc.). In this case the 'telegraph equations' are joined into a system with correspond- 
ing boundary value conditions. All the existing methods of calculating the currents 
in these systems are based on finite-dimensional approximation of the continuous long 
lines (finite-elements methods, etc.), and so there is a loss of precision in this approxima- 
tion. These methods usually need a prior manual analysis of the network topology. For 
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some particular kinds of problem (when interaction between the processes at the differ- 
ent branches is realized only in a form of Kirchhoff 's law at the nodes where branches 
are connected), we obtain a new method based on the approach [2] for the partial dif- 
ferential equations on graphs. We derive explicit formulas for the Fourier transforms 
of the currents and voltages in the branches of the branching network. These formulas 
give an exact solution of the system of 'telegrapher's equations'; the corresponding al- 
gorithm is suitable for computer applications. The input data are provided directly by 
the network topology, without the necessity for prior manual analysis to obtain addi- 
tional equations for boundary- value description. The availability of the formulas for the 
currents is used to derive the explicit formulas for currents amplitude maximum being 
achieved over all the external influences (inputs) satisfying the constraints for input 
energy or input value at every point. This maximum may be used as a characteristic 
of network robustness to excess voltage. 

Let us consider a system of n long lines with lengths h, . . . , l n that are connected to 
a network (graph with branches and nodes). The equation for the current and voltage 
in each branch with number k = 1, . . . , n is [1] 

m£{y,u) = -(R k + iujI k )I k (y,uj) + e k (y,u) 

' J _ (!) 

-5${y,u) = -(Y k + tu;C k )U k (y,uj) + 5(y - y k )I k (uj) 

In eqn.l, I k (y,uj), U k (y,uj) are the Fourier transforms of current and voltage in the 
kth branch, y G [0, l k ] is a longitudinal co-ordinate (we have chosen the orientation 
for every branch), u> is an angular frequency, to G R. The constants R k , I k are the 
series resistance and inductance per unit length. The constants Y k , C k are the parallel 
conductance and capacitance per unit length. The action of the external electric field 
is characterized by the functions e k (y,uj); I k (u>) are the external currents flowing into 
the given points y k G [0, I k \; and 5 is a delta function. 

We assume that interaction between the processes at the different branches is re- 
alized only at the nodes where branches a re connected, and so all the eqns. 1 for 
k = 1, . . . , n are joined into the system with supplementary boundary- value conditions 
in the form of Kirchhoff 's current and voltage laws at joining nodes. Ohm's law holds 
at the terminal nodes that belong to only one branch, so that I k (0, uj) = Y k U(0,uj) or 
I k (I k: uj) = Y k U k (I k ,uj), where the constants Y k are the terminating admittances. 



2 



We suppose that 

(i) Nonzero Rk or exist for some k. 

(ii) \ek(y,uj)\ < c for some constant c > for every k, y, u>. 

2 Formulas for currents and voltages 

Our method is taken from Reference 2. The main idea is to describe the current and 
voltage distribution by the complex 2n-vector 



defined for every x G [0, l],u> € M.. This vector is convenient because we can write the 
boundary value problem for it: 

dz 

— — (x, to) = A(x, lo)z(x, uj) + E(x, u), 



Here the external influence is 

E(x,lo) = [hei(h,uj),5(lix - yi)Ji(cj), . . . , e n (l n x,u), S(l n x - y n )I n (uj)} T (4) 

In eqn.3, A(uj) = diag [A\ (u) , . . . ,A n (u)] is a complex 2n x 2n-matrix, where 2x2 
matrices 



This means that the matrices are placed on the main diagonal of A and all the 
remaining elements of A are zero. 

The In x 2n-matrices Bq and B\ contain complete information about the topology 
of our network; we shall give the algorithm for their construction below: 

The explicit final formula for the Fourier transform of the current and voltage is 



z(x,lo) = [Ui{hx,Lo),Ii(hx,uj), U n (l n X,Uj),I n (l n X,Uj)} 



T 



(2) 



B o z(0,u) +B 1 z{l,u) = 0. 



(3) 





(5) 



together with eqn.2. 
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In eqn.5, G(x, p, 00) is the Green's function for the problem of eqn.3 defined [3] for 
x G [0, 1], p G [0, 1], w 6 K as a complex 2n x 2n matrix 



G(x,p,u) 



-$(x,w)[B + 5i$(l,a;)] _1 .Bi$(l - p,uj) for < x < p 
$(x,u)[B + 5i$(l,u;)]- 1 5o$(-p,u;) for p < x < 1 



(6) 



Note that the 2/cth component of eqn.5 is undefined at the point x = l^kVik-, an d 
has a jump discontinuity at this point for k = 1, . . . , n (see eqn.l). 

The matrix <&{x,u}) = exp[A(u)x\ is the fundamental matrix of the first eqn.3: 
&(x,u>) = diag [3>i(cc, 00), ... , <3? ra (x, u)], where 



cosh(7 fc x) -g- sinh(7 fc x) 
-gsinh(7 fe x) cosh(7 fc x) 

= h(Rk + ivL k ), ( k = h(Y k + iivC k ), Jk = V^kCk (we take the principal value for 
the root and assume that sinh 0/0 = 1). 

The topology of the network is defined by the ordered set 
{s m ,p m ,r m , . . . ,i Pm ), (ji, ■ ■ ■ ,jr m )}m=i- Here 1 is tne number of nodes, p m 
is the number of branches with numbers i±, . . . ,i Pm existing from the mth nodes, r m 
is the number of the branches with numbers ji, ■ ■ ■ ,j Pm entering into the mth nodes, 

Sm = Pm f 'rn- 

We construct the matrices Bq and B\ by joining all the strings of s m x n-matrices 

Qm = Qm(i,j), Prn = P(i,j),m = 1, . . . , q, 1 = 1, . . . , S m , j = 1, . . . , TV. 





Pi 




Qi 


B = 




Bi = 










Qq 



where Q m and P m are constructed by the following algorithm. 

(i) Let s m = 1. If p m = 1, then P m (l, 2i ± - 1) = Y k and P m (l, 2ii) = -1; if r m = 1 
then Q m (l, 2ji — 1) = Yk and Q m (l, 2j\) = —1. All the remaining components of 
P rn and Q m must be zero. 

(ii) Let s m > 1- If r m > 0, then Q m {Pm + d, 2j d -l) = -Q m (p m + d, 2j d+1 - 1) = 1 for 
d = 1,... ,r m - 1, and Q m (s m ,2ji - 1) = Q m (s m ,2j 2 ) = ■ ■ ■ = Qm{s m ,2j rm ) = 1. 
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If p m > 0, then P{s,2i s - 1) = -P m (s,2i 8+1 - 1) = 1 for s = 1, . . . ,p m - 1 
and P m (s m ,2i 2 ) = ••• = im(sm, 2i p0 o) = 1- If r m > 1 and p m > 1, then 
Pm(Pm, 2i P m — 1) = —Qm{Pm, 2ji — 1) = 1. All the remaining components of P m 
an Q m must be zero. 

Let us discuss the correctness of our formulas. The existence of nonzero Rk or Y k is 
sufficient for matrix [Bq + Bi$>(l,u)] to be invertible, and for the problem of eqn.3 to 
be well posed [3]. Let us show this. Let zq G C™ and [Bq + Bi&(1,uj)]zq = 0. There is 
an energy dissipation in eqn.l for k with nonzero Rj, or Y^, and so the corresponding 
evolution system is stable in the temporal domain. Thus zq = 0, and the problem in 
eqn.3 is well posed. 

Integral eqn.5 exists because of supposition (ii). 

3 Maximum of current amplitude 

Having obtained explicit formulas for currents, we can derive the explicit formulas for 
the maximum current amplitude at a branching network. This problem is important for 
applications, because we can use this value for the characterization of the robustness to 
excess voltage or electromagnetic impulse. We think that this problem has no explicit 
solution other than our explicit formulas for currents. 

In eqn.l, the external influences are characteristised for every to 6 K by the ordered 
set F(u) = {e k (.,uj),I k (uj),y k }l =l . 

Fix a k > 0,b k > 0, < < < l k for k = 1, . . . , n. Let us introduce the sets 
Fp(oj) = {Tiuj)} being the classes of admissible external influences: 

(i) For 1 < p < +oo, class F p {ui) is the set of such T(u>) = {] y (., u), T\\ (oj), f y }n =00 



that 




< a k , \I k {uj)\<b k , y k G [a k ,(3 k ], k = l,...,n. 



(ii) Class Foofu;) is the set of such F{oj) = {] y (., u), X\\ (u>), t||}n =00 that 



sup \e k (y,u)\ < a k , \I k (u)\ < b, 
ye[o,i k ] 



> k , y k e [a k ,/3 k ], k = l,...,n. 
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Our aim is to obtain the value 



sup \I k {y,u)\ 



V 



for every k = 1, . . . , n, y € [0, 1 < p < +oo. 
The result is that, for p £ (1, +oo], we have 



sup = V 

^M^M m=l 



afc 



/ |G2fe,2m-l(2/Afc,P,w)| 9 dp 

vo 



1/9 ^1 

+ sup |G2fe, 2m (y//fe,p,w)| > (7) 



where g = p(p — 1) 1 for every p G (1, oo), g = 1 for p = +oo. For p = 1, we have 

sup = X] 

< Ofc sup |G 2 fc,2m-l(yAfc,P,w)| +6fc sup |G 2 fc,2m(yAfc,P,w)| ? (8) 

[ pe[o,i] pe[a k /ik,Pk/ik] ) 

Here, Gk >m are components of matrix G defined by eqn.6. 

Eqns.7 and 8 are very suitable for computer application. Analogous results can be 
obtained for 

sup \U k (y,u)\ 
F(u>)eT (w) 

V 

(the corresponding changes in eqns.7 and 8 are obvious). 

To prove eqns.7 and 8, we have to remark that the right-hand part of eqn.7 or eqn.8 
is the functional norm of the 2/cth string of G(y/l k , ■, w) being presented as an element 

of a Banach space [4], which is dual to such space of function e k (-, u>) that T{(S)eT (u) 

V 

defines a ball there. 
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